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Abstract 



O 

CN ' Background: Inelastic neutrino-nucleus scattering through the weak neutral-current plays im- 

b 

portant role in stellar environment where transport of neutrinos determine the rate of cooling. Since 

there are no direct experimental data on neutral-current neutrino-nucleus cross sections available, 
CO 

only the modeling of these reactions provides the relevant input for supernova simulations. Pur- 
pose: To establish fully self-consistent framework for neutral-current neutrino-nucleus reactions 
based on relativistic nuclear energy density functional. Methods: Neutrino-nucleus cross sec- 
tions are calculated using weak Hamiltonian and nuclear properties of initial and excited states are 
obtained with relativistic Hartree-Bogoliubov model and relativistic quasiparticle random phase 



£> ■ approximation that is extended to include pion contributions for unnatural parity transitions. Re- 

suits: Inelastic neutral-current neutrino-nucleus cross sections for 12 C, 16 0, 56 Fe, 56 Ni, and even 

(N 

isotopes 92_100 Mo as well as respective cross sections averaged over distribution of supernova neu- 

co ' 

^^ ■ trinos. Conclusions: The present study provides insight into neutrino-nucleus scattering cross 



sections in the neutral channel, their theoretical uncertainty in view of recently developed micro- 
scopic models, and paves the way for systematic self-consistent large-scale calculations involving 



H ■ open-shell target nuclei. 
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I. INTRODUCTION 



Nuclear weak interaction processes play an important role in the evolution of supernova 
collapse, e.g., electron capture, beta-decay, and neutrino- nucleus reactions (l|. Neutrino- 
nucleus scattering through the weak neutral current could provide contributions of relevance 
in stellar environment where transport of neutrinos determine the rate of cooling |l|-|3j. 
Recently, inelastic neutrino-nucleus scattering has been introduced in supernova simulations 
as novel mode of energy exchange between neutrinos and matter 4|-|6J. Although this process 
has no large effect on collapse trajectories, it has significant contribution to increasing the 
neutrino opacities, and it strongly reduces the high-energy tail of the neutrino spectrum 
emitted in the neutrino burst at shock breakout [5|. To date, only a single microscopic 
framework for the neutral-current neutrino-nucleus scattering, based on hybrid model 7|, 
has been included in supernova simulations (l||. Since the calculations of weak interaction 
processes in various theoretical models can result in differences in the reaction rates and 
cross sections, sometimes larger than order of magnitude |8J-lll). providing insight into the 
neutral-current neutrino-nucleus cross sections from independent models is paramount for 
assessing the impact of the uncertainties in nuclear structure models on the outcomes of 
supernova simulations. 

Modeling neutrino-induced reactions is also important in view of studies on modern detec- 
tors based on neutrino scattering on hadrons and nuclei. The ongoing and planned neutrino 



detector facilities involve variety of target materia 



fives. These include MOON 
OPERA 



12|], MiniBooNE_p|, NEMO 



s, induced reactions and scientific objec- 
MINOS (l5|, SNO+ (l6|, 



iica 

Q. 



17fl, LVD (Large Volume Detector) 18], ORLaND experiment at the Spallation 



Neutron Source (SNS) 19J . NOvA neutrino experiment [20[, etc. In addition to charged- 
current neutrino-induced processes in detectors, another possible reaction channel includes 
neutral-current neutrino scattering on hadrons and nuclei, resulting e.g., in small showers 
of secondary gamma radiation and creation of electron-positron pairs. Although the cross 
sections in the neutral channel are smaller than in the case of charge-exchange reactions, 
understanding complete detector response necessitates consistent microscopic insight into 
all relevant processes involved. 

Over the past years, several theoretical frameworks have been developed to provide de- 
scription of the inelastic neutrino-nucleus scattering in the neutral channel. Due to consider- 



able progress in the shell model Hamiltonians, a number of neutrino-induced reactions have 



been described, also including various decay channels 



2lN26|. Random phase approximation 



(RPA) based on Landau-Migdal force has been employed in calculations of neutrino-induced 



27 



28|. The hy- 



reaction rates for r-process nuclei, including those in the neutral channel 

brid model combines the shell model for allowed transitions, with the RPA to account for 

the forbidden transitions, allowing systematic calculations for a large number of target nu- 



clei 



IT. 



29 



30J. The Hartree-Fock + RPA based fully on 



Sk yrm e functional has been 



31 



32] . In view of developing 



employed in studies of the cross sections for 12 C, 16 0, and 208 Pb 
neutrino detectors and astrophysical role of neutrino-nucleus reactions, the relevant nuclear 
matrix elements have recently been revisited in Ref. 33|, and employed in studies involving 
target nuclei 40 Ar, 56 Fe, 92 " 1()0 Mo and 128 < 130 Te, based on quasiparticle RPA (QRPA) |3J-|37_|. 



In another recently developed framework based on Hartree-Fock-Bogoliubov (HFB) model 
and QRPA, Brueckner G matrix is employed for two-body interaction by solving the Bethe- 
Salpeter equation based on the Bonn CD potential, and pairing correlations have also been 



taken into account 



38H40J. Supernova neutrino- Fe cross sections in the neutral channel 



have also been explored in the local density approximation taking into account the Pauli 



blocking and Fermi motion effects 



4l| 



Whereas for the charged- current neutrino-nucleus reactions some experimental data are 
available for 12 C and 56 Fe 42 



in the case of neutral- current inelastic neutrino-nucleus 
scattering there are no experimental data available, except for the ground state transition 
to the 15.11 MeV state (T=l) in 12 C |46l . |47| . The indirect experimental insight into the 
inelastic neutrino-nucleus cross sections can be obtained from inelastic electron scattering. 



As shown in Ref. 



48|, magnetic dipole strength distributions for several iron group nuclei 



are dominated by isovector Gamow- Teller transitions that can be translated into inelastic 
neutral-current neutrino-nucleus cross sections. Due to lack of direct experimental data on 
these cross sections, modeling by various approaches is crucial. In this way one can provide 
not only the relevant input for supernova simulations, but also the insight into theoretical 
uncertainties in description of the cross sections in the neutral channel. 

In this paper we introduce the framework for the neutral-current neutrino-nucleus inelas- 
tic scattering based on relativistic nuclear energy density functional. Within this framework 
the nuclear ground state and various excitations induced in nuclei by the incoming neutri- 
nos are described in a fully self-consistent approach, i.e., universal effective interaction is 



employed without introducing any adjustments to the specific properties of target nuclei or 



neutrino energies involved. Although the energy density functional has a 
ployed in the non-relativistic framework using Skyrme parameterizations 31 



ready been em- 



32fl , the present 



study provides the first self-consistent framework to describe neutrino-induced reactions in 
the neutral channel, involving open-shell target nuclei that necessitate explicit implementa- 



tion of the pairing correlations. The relativistic nuclear energy density functiona 



successfully employed in studies o 
decay rates of r-process nuclei 



giant resonances and exotic modes of excitation 



has been 



49 



-54 



55] , muon capture |56j and stellar electron capture rates [11 1 



and in constraining the neutron skin in nuclei 



57 



In Refs. 



um 



the relativistic proton- 



neutron QRPA has been employed in modeling charged-current neutrino-nucleus reactions. 
In the present analysis of neutrino-induced reactions in the neutral channel, the model ne- 
cessitates further development of the relativistic QRPA outlined in Ref. 60| in order to allow 



taking into account both natural and unnatural parity excitations in the neutral channel. 

The paper is organized as follows. In Sec. II we introduce the basic formalism for the 
neutrino-nucleus cross sections in the neutral channel based on weak Hamiltonian and rela- 
tivistic nuclear energy density functional. The respective cross sections have been explored 
in detail for a set of target nuclei in Sec. III. The conclusions of the present work are 
summarized in Sec. IV. 



II. THEORETICAL BACKGROUND 



The weak process to be considered is inelastic neutral-current neutrino-nucleus reaction, 



v e + z X N ->• v, 



~z A 



n- 



(1) 



where the incoming electron neutrino (u e ) scatters on target nucleus X(Z, N) which absorbs 
part of the neutrino energy. The interaction between the neutrino and nucleus is described 
by weak Hamiltonian, while the properties of initial and final states of target nucleus are 
described by effective nuclear interaction, in this particular case formulated using relativistic 
energy density functional. The formalism leading to the expression for the cross section is 



given in Refs. 



61 



62| . The general expression for neutrino-nucleus differential cross section 



is derived in terms of relevant multipoles of the nuclear weak currents, 
da u \ G\(? At: 



r F < 



dVL J 2ir 2 2Ji + l 



I J>0 

-2q(u + f3)Re(J f \\£ J \\J i )(J f \\Mj\\J i y} 

+ ^{(l-(r>-g)(/3-4))[|(J/||f/ AG ||J J )| 2 + |(J / ||f/ L ||J J )| 2 " 
j>i 

+2q-( I >-/3) J Re(J / ||f/ /AG ||J J )(J / ||f/ L ||J J )*}|, (2) 

where Gp is Fermi constant for the weak interaction and e denotes the energy of outgoing 
neutrino. The momentum transfer q = v — k is defined as the difference between the 
incoming (u) and outgoing (fe) neutrino momenta, q and u denote the corresponding unit 
vectors, and (3 = k/e. The transition matrix elements between the nuclear initial and 
final states include transition operators of various multipoles: charge Aij, longitudinal Cj, 
transverse electric 7j L , and transverse magnetic 7} M AG multipole operators, expressed in 
terms of spherical Bessel functions, spherical harmonics, and vector spherical harmonics 



611, 



62|. Complete calculation of inelastic neutrino- nucleus scattering necessitates inclusion of a 
number of multipoles J . Although higher-order multipoles have rather small contributions 
at low incoming neutrino energies, these can not be neglected at energies about tens of 



MeV 



31 



32 



40|. In the present study, multipoles up to J = 5 contributing to the cross 
section in Eq. ([2]) will be included in calculations. In the specific case of the neutrino-nucleus 
scattering in the neutral channel, the transition operators listed above include the following 



form factors 



33. 



63| 



vector form factors F?, fi v = F^ - 2MF^ 



^^-^{^mkvf) (3) 

^> tf) = I (1 -W M (l + IS 4 ¥ )" (4) 

^ (f) ^v^^ l+ _A_y (5) 



axial vector form factor 



F ^\ f) ^\ Fm {, + J -^— f] (6) 



pseudoscalar form factor 

9 i 9 

qz -\. m %. 



fP-% 2 ) = "T^J" ' ■ < 7 > 



The indices n,p denote the respective form factors for neutrons and protons, 6w denotes 
the Weinberg angle, sin 2 8 W =0.2325, and the static values are F A (0) = -1.2617, /i y(n) (0) = 
-1.463, and //^^(O) =1.054. In the present analysis the strange quark content in the form 
factors has been neglected. By employing the full operator structures in the transition 
matrix elements, the inelastic neutrino- nucleus cross section is evaluated using Eq. (J2J), with 
an additional quenching factor included in the free-nucleon axial- vector coupling constant, 
resulting in -Fa(O) = 1.0. This quenching corresponds to additional factor 0.8 in Fa (EJ). 
The value of the quenching depends on the effective interactions and the model space under 
consideration, e.g., in the hybrid model quenching factor of 0.74 has been used |7J|, while 0.8 
represents reasonable value for the framework employed in the present study [9|. 

The transition matrix elements between the initial and final states in Eq. (J2J) are deter- 
mined in a fully self-consistent framework based on relativistic nuclear energy density func- 
tional 64j, |65j. Therein the nuclear ground state is described with the Relativistic Hartree- 



Bogoliubov (RHB) model, and excited states are calculated using the relativistic quasipar 



tide random phase approximation (RQRPA) 5l|, |60|. The application of relativistic nuclear 



energy density functional is realized in terms of the self-consistent mean field theory for nu- 
cleons and minimal set of meson fields; isoscalar scalar a-meson (J n = + ,T = 0), isoscalar 
vector w-meson (J w = 1~,T = 0) and the isovector vector p-meson (J 77 = 1~,T = 1), 
supplemented with the electromagnetic field. The meson-nucleon interaction is included 
with a minimal set of the interaction terms, where the vertex functionals include explicit 
dependence on the vector density. The details of the RHB model based on this class of 



effective density- dependent interactions are given in Ref. 66J). For the model parameters 
that determine the density-dependent coupling strength and the meson masses we employ 
the values of the DD-ME2 parameterization, obtained by simultaneous adjustment of the 
effective interaction to the binding energies, charge radii, differences between radii of neu- 
tron and proton density distributions for 12 spherical nuclei and nuclear matter properties 

6 



at saturation density 67[. The pairing correlations in open shell nuclei are described by the 



finite range Gogny interaction, with parameterization D1S 



68|. 



The RQRPA is formulated in the canonical single-nucleon basis of the RHB model, and 
the residual interaction is derived from the same nuclear energy density functional as in the 



RHB model 



60 



It includes not only configurations composed from two-quasiparticle 
states of positive energy, but also pair-configurations formed from the fully or partially 
occupied states of positive energy and empty negative-energy states from the Dirac sea. 
In the implementation in modeling the weak interaction processes, the major advantage of 
the RHB+RQRPA model is that it is fully consistent in view of the effective interactions 
employed. In the particle-hole (ph) and pairing (pp) channels, the same interactions are used 
in the RHB equations that determine the canonical quasiparticle basis, and in the matrix 
equations of the RQRPA. In this way, one can employ the same nuclear energy density 
functional in description of the weak processes throughout the nuclide map without any 
additional adjustments of the model parameters. 

In the present study we further extend the RQRPA framework outlined in Ref. |6C 
by including the pion contributions in order to account both for the natural, (— 1) J = ir, 
and unnatural, (— 1) J+1 = n, parity excitations that take part in inelastic neutrino-nucleus 



scattering. As shown in Refs. 



38 



40[, unnatural parity excitations play an important role 
in the overall neutrino-nucleus cross sections in the neutral channel, and should be included 
in a framework aiming to provide consistent and reliable results. Excitations of unnatural 
parity states necessitate the inclusion of the residual interaction term generated by the it- 
meson (J 77 = 0~,T = 1) exchange. At the Hartree level, i.e. in the RHB, the pion does not 
contribute because it carries unnatural parity and the corresponding mean field breaks parity. 
The pion major effect comes from the second and higher order diagrams in the correlated 
two-pion exchange. The quantum hadrodynamics model (QHD II) included in addition to 
(a,u;,p) meson fields also a pseudoscalar pion field. However, as pointed out in the RRPA 
study in Ref. [70), the pseudoscalar pion couples too strongly, resulting in total disruption of 
the ordering of the lowest excited states. The RRPA analysis showed that implementation 
of pseudovector pion-nucleon coupling improves the spectrum in comparison to experiment, 
especially for the pion-like states J w = 0~,2~,4~. Therefore, in the present study we use 



the pion-nucleon interaction with pseudovector coupling, given in the Lagrangian density as 

4 PV) = — ^sT^tttV (8) 



and the propagator for the residual two-body interaction reads, 

D^\q) = - 2 1 - . (9) 

The standard value for the pseudovector pion-nucleon coupling is / 2 /47r = 0.08, while the 
measured pion mass amounts m n = 138 MeV. Since the one-boson-exchange interaction 
with pseudovector coupling (jSJ) contains a__contact term, one accounts an additional term 
for the 5-force to remove its contribution 



The two-body matrix elements of the one- 
pion exchange interaction and the 5- force in pseudovector coupling are calculated in the 



momentum space representation according to detailed formalism given in Ref. [7l| . When 
calculating the neutrino-nucleus cross sections in Eq. ([2]), for each transition operator Oj 
the matrix elements between the ground state and the final state of target nucleus are 
expressed in terms of single-particle matrix elements between quasiparticle canonical states, 
the corresponding occupation factors v^^u^ and forward- and backward-going amplitudes 



60 1 



X, Y, obtained by diagonalization of the RQRPA matrix 

(JfWOjU) = Y, {^MIOjIIaO + (- 1 ) JM " V+J ^ (A\6j\\v) } 

x(tW+ (-l)%v)- ( 10 ) 

All relevant transitions between the |0 + ) ground state and \J* ) final states are taken into 
account in the following calculations. 

III. RESULTS AND DISCUSSION 

We have employed the framework introduced in Sec. II in modeling the neutral-current 
neutrino- nucleus scattering for a set of target nuclei of interest for neutrino detector response 
and understanding the role of neutrinos in supernova evolution. In particular, the cross 
sections have been calculated as a function of the incoming neutrino energies for 12 C, 40 Ar, 
56 Fe, 56 Ni, and 92_100 Mo isotopes. The nuclear matrix elements are obtained using the 
energy density functional with DD-ME2 parameterization 67J, supplemented by the Gogny 



force D1S to account for the pairing correlations in open shell nuclei 



68]. The overall cross 



section Eq. ([2]) includes summation over transitions to all possible final states characterized 
by multipoles up to J = 5 with both positive and negative parity. 

In Fig. [IJthe calculated electron neutrino-nucleus cross sections are showed for the inelastic 
scattering 12 C(z/ e , z/g) 12 C* for the range of neutrino energies ^=0-100 MeV. Complete set 
of multipole states J 71 = 0^ — 5^ is taken into account in the overall cross section. Various 
contributions from the most relevant multipole states J 77 = ± — 3 1 * 1 to the cross sections are 
also displayed separately. For comparison, Fig. [1] also shows recent results for the overall 



38]. A 



s one can 



cross sections for 12 C target, based on QRPA with Bonn CD potential 
observe in this figure, at low neutrino energies the overall cross sections are dominated by 
1 + transitions. However, as the energy increases to 100 MeV, the role of other multipoles 
becomes important, in particular those of 1~, 2~ and smaller contribution from 2 + states. 
In the specific case of neutrino energy E u =50 MeV, the present results are at variance with 
Ref. |3jJ, where multipole contribution from J n = 1~ dominates over 1 + . On the other 
hand, the multipole composition of the cross sections is in qualitative agreement with recent 
study based on the QRPA 38] . The total cross sections from the present study appear 



38|. We turn to 



systematically larger than the respective values obtained using the QRPA 
this discrepancy later in discussion of the cross sections. 

In Fig. [2] the cross sections are shown for the scattering process 40 Ar(^ e , i/) 40 Ar*. In 
comparison to 12 C, the interplay between various multipoles becomes more involved. Al- 
though at low-energies 1 + transitions dominate, at E v above ^40 MeV, 1~ and 2~ multipoles 
have the largest contributions. Neutrino-induced reactions with 40 Ar have been studied in 



details in recent work based on QRPA 



40] . in view of their relevance for detecting core 



collapsing supernovae neutrinos. For comparison, the respective QRPA results from (40j are 
also shown in Fig. |2j The total cross sections from the present analysis appear up to an 



order of magnitude larger than the QRPA 



40] ones. Even though a variety of advanced 



theoretical frameworks have been developed over the past years, one can observe consider- 
able theoretical uncertainty inherent in the modeling of the neutral-current neutrino-nucleus 
cross sections. These uncertainties originate to a large extent to differences in single-particle 



spectra and respective transitions induced by incoming neutrinos. In Ref. J32J it has been 
shown that even within the same model, Hartree-Fock + RPA based on Skyrme functional, 
only small adjustment of single-particle parameters resulted in 30% increase of the overall 
neutral-current neutrino-nucleus cross sections. Therefore, it is not surprizing that imple- 



mentation of various models using independent effective interactions could result in large 
differences. Consequently, it is important to provide reasonable quantitative estimate of 
theoretical uncertainty in the cross sections and to critically assess its effect in modeling 
supernova evolution and neutrino detector response. In recent analysis of charged- current 
neutrino- nucleus cross sections |9|, by employing variety of microscopic models and effective 
interactions, it has been shown that one can provide reasonable estimate of the theoretical 
uncertainty in modeling weak interaction processes. 

Figures [3] and H] show the neutral-current neutrino- 56 Fe and - 56 Ni cross sections, re- 
spectively. Although the multipole composition of the cross sections appear in qualitative 
agreement, some smaller differences can be noted due to differences in neutron and proton 
numbers and respective excitation spectra. However, one can conclude as general property 
that J = 1 states are the most dominant, at lower energies J w = I + dominates while at 
E v <J 65 MeV transitions J w = I - have the major contribution. In addition, at higher ener- 
gies J n = 2~ state also competes with J n = 1~ for dominance. At high-end neutrino energy 
ps 100 MeV other multipole transitions also contribute to the overall cross sections, e.g., 
jit _ 2± 5 3±^ ^he calculated cross sections for the scattering process 56 Fe(^ e , z/) 56 Fe* are ex- 
plored in detail in comparison with the hybrid model 30| and QRPA based framework |34j. 



In Table [J the RQRPA cross sections are given for a selection of neutrino-energies up to 100 
MeV. For comparison, the results are shown both with and without quenching in Fa (EJ). It 
is interesting to observe that at low neutrino energies, where the cross sections are rather 
sensitive on the fine details of the transition spectra, the RQRPA (with quenching) and 
hybrid model results are in excellent agreement. Although the cross sections from the three 
models appear in the overall qualitative agreement, in the energy region of relevance for the 
supernova neutrino processes (ps20-40 MeV) the RQRPA cross sections are up to a factor 
psl.5 (2.0) larger than the hybrid model and QRPA results, respectively. It is interesting 
to note that very recent QRPA study resulted in averaged cross sections for 56 Fe roughly a 
factor of two larger than for the hybrid model 37 1. 



The scattering cross sections in the neutral channel have also been explored for a set 
of Mo isotopes, that recently became interesting due to on-g oing and future applications 
of molybdenum in terrestrial neutrino detectors, MOON 12j and NEMO |14J . related to 



neutrino studies and search for the events of neutrinoless double beta decay. In the present 
analysis the cross sections have been explored for the most abundant even molybdenum iso- 

10 



E v [MeV] w/o quench, w quench. Hybrid [30] QRPA [34] 



10 


2.91(-1) 


1.87(-1) 


1.91(-1) 


1.01(4-0) 


15 


4.30(+0) 


2.77(+0) 


2.19(4-0) 


2.85(4-0) 


20 


1.51(4-1) 


9.78(+0) 


6.90(4-0) 


5.79(4-0) 


25 


3.42(+l) 


2.22(+l) 


1.51(4-1) 


1.06(4-1) 


30 


6.26(+l) 


4.08(+l) 


2.85(4-1) 


1.87(4-1) 


35 


1.04(+2) 


6.79(+l) 


4.89(4-1) 


3.24(4-1) 


40 


1.57(+2) 


1.05(+2) 


7.86(4-1) 


5.51(4-1) 


45 


2.32(+2) 


1.54(+2) 


1.19(4-2) 


9.05(4-1) 


50 


3.24(+2) 


2.16(+2) 


1.72(+2) 


1.43(+2) 


55 


4.38(+2) 


2.94(+2) 


2.39(+2) 


2.15(4-2) 


60 


5.76(+2) 


3.89(+2) 


3.20(+2) 


3.09(+2) 


65 


7.40(+2) 


5.01(+2) 


4.15(4-2) 


4.26(+2) 


70 


9.29(+2) 


6.33(+2) 


5.25(+2) 


5.63(+2) 


75 


1.15(+3) 


7.85(+2) 


6.50(+2) 


7.17(+2) 


80 


1.40(+3) 


9.59(+2) 


7.89(4-2) 


8.82(+2) 


85 


1.68(+3) 


1.16(+3) 


9.42(+2) 


1.05(4-3) 


90 


2.00(+3) 


1.38(4-3) 


l.ll(+3) 


1.22(4-3) 


95 


2.36(+3) 


1.63(4-3) 


1.29(4-3) 


1.38(4-3) 


100 


2.76(+3) 


1.92(4-3) 


1.49(4-3) 


1.52(4-3) 



TABLE I. The cross sections for 56 Fe(z^ e , fg) 56 Fe* process, given in units of 10 42 cm 2 . The results 
of the present analysis without (second column) and with th e q uenching factor (0.8) in Fa (third 



column) are com par ed with the results of the hybrid model 



model from Rcf. 



30fl (forth column) and QRPA based 



3J (fifth column). 



topes, 92 Mo, 94 Mo, 96 Mo, 98 Mo and 100 Mo. The major contribution in natural molybdenum 
comes from 98 Mo, amounting 24.13%. Table [Ln shows the neutral-current neutrino-nucleus 
cross sections for even isotopes 92_100 Mo in the range of incoming neutrino energy £^=10- 
100 MeV. By inspecting the numbers, one can observe rather small but systematic increase 
in the cross section values for all neutrino energies. As can be expected, 100 Mo has the 



11 



E„ [MeV] 92 Mo 



91 



Mo 



<)(, 



Mo 



<)X 



Mo 



100 



Mo 



10 

15 
20 
25 
30 
35 
40 
45 
50 
55 
60 
65 
70 
75 
80 
85 
90 
95 
100 



6.32(-2) 3.58(-l) 5.25(-l) 6.67(-l) 7.82(-l] 



2.33(+0 


) 3.24(+0 


3.67(40 


4.01(40 


4.26( 


9.68(+0 


) 1.14(+1 


1.21(41 


1.26(41 


1.29( 


2.47(+l 


) 2.76(+l 


2.86(41 


2.95(41 


3.00( 


4.98(+l 


) 5.42(+l 


5.59(41 


5.73(41 


5.84( 


8.95(+l 


) 9.61(+1 


9.88(41 


1.01(42 


1.03( 


1.46(+2 


) 1.55(+2 


1.59(42 


1.63(42 


1.67( 


2.23(+2 


) 2.36(+2 


2.42(42 


2.48(42 


2.54( 


3.22(+2 


) 3.39(+2 


3.47(42 


3.56(42 


3.64( 


4.44(+2 


) 4.66(+2 


4.78(42 


4.89(42 


5.00( 


5.89(+2 


) 6.16(+2 


6.32(42 


6.47(42 


6.61( 


7.58(+2 


) 7.91(42 


8.10(42 


8.30(42 


8.47( 


9.49(+2 


) 9.88(42 


1.01(43 


1.03(43 


1.06( 


1.16(+3 


) 1.21(43 


1.24(43 


1.26(43 


1.29( 


1.40(+3 


) 1.45(43 


1.48(43 


1.51(43 


1.54( 


1.65(+3 


) 1.71(43 


1.75(43 


1.79(43 


1.82( 


1.92(+3 


) 1.99(43 


2.03(43 


2.07(43 


2.11( 


2.22(+3 


) 2.29(43 


2.34(43 


2.38(43 


2.42( 


2.52(+3 


) 2.61(43 


2.66(43 


2.71(43 


2.75( 



-o) 

-1) 
-1) 
-1) 

-2) 
-2) 
-2) 
-2) 
-2) 
-2) 
-2) 
-3) 
-3) 
-3) 
-3) 
-3) 
-3) 
-3) 



TABLE II. The total neutral-current neutrino-nucleus cross sections for even isotopes 92 100 Mo, 
given in units of 10 2 cm 2 . 



largest cross section of all shown so far simply by virtue of large number of active nucleons 
contributing to the collective nuclear response in the scattering process. The respective 
cross sections for 98 Mo and their multipole composition are displayed in Fig. For neutrino 
energies below 10 MeV both 0~ and 1 + states have relevant contributions. At E u « 45 MeV 
one can observe the intersection between the main components in the cross sections: 1 + at 
lower energies and 1 _ , 2~ which dominate at higher energies. When comparing the overall 



cross sections to those of other recent studies based on QRPA 



371 140| . the present results 
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are systematically larger, but within an order of magnitude. 

In order to explore theoretical uncertainties in modeling neutral-current neutrino-nucleus 
reactions in more details, in Figs.[6]and[7]partial multipole contributions to the cross sections 
for 96 Mo target are shown at incoming electron neutrino energies £^=20 and 100 MeV, 
respectively. The results of the present study (RQRPA) are shown in comparison with the 



cross sections recently obtained using QRPA (Balasi et al., Ref. |37fl). Although the total 
RQRPA cross sections are somewhat larger than those of QRPA, one can observe to a large 
extent excellent qualitative agreement between the two models based on rather different 
backgrounds. At low neutrino energy (Fig. [6]) in both cases largely dominant excitation 
channel is 1 + . The distribution over various multipoles appears rather involved at E Ue =100 
MeV. The main contribution is obtained for 1~ transitions, but other multipoles also show 
considerable effects, ranked in the order of importance as follows: l - , 2 + , 2~, 1 + , 3 + , 3~, 
etc. The models based on RQRPA and QRPA result in excellent agreement in relative 
contributions of various multipoles, except for the anomaly for the QRPA 1 + channel. 

An important application of microscopic models of neutrino-nucleus reactions is descrip- 
tion of the cross sections for stellar neutrinos of relevance for the neutrino detectors that 
could provide better insight into fascinating events in the universe that produce neutrinos. 
The calculated cross sections given as functions of the incoming neutrino energy can be aver- 
aged over supernova neutrino flux, that is usually described by the Fermi-Dirac distribution, 

1 E 2 

IXEu) = T*exp[(E v /T)-a] + l " (U) 

Especially interesting is modeling the reaction rates of neutrinos scattering on nuclei that 
can be used as targets for the supernova neutrino detectors, e.g., 40 Ar, 56 Fe, 56 Ni, Mo 
isotopes, etc. In this way, one can predict expected number of events in detector that 
originate from specific stellar environment which determines the production of low-energy 
neutrinos. In this work we calculate the neutral-current neutrino-nucleus cross sections 
averaged over the supernova neutrino flux in the range of temperatures T v = 2 — 10 MeV, 
and for the chemical potential a=0. Figure |S] shows the respective flux-averaged cross 
sections for a set of target nuclei, 12 C, 40 Ar, 56 Fe, 56 Ni and 98 Mo. As the temperature 
increases, neutrinos with higher energies have larger contributions in the averaged cross- 
sections. The reason is two-fold, i) the Fermi-Dirac distribution shifts toward higher energies 
with increased temperature, and ii) the neutrino-nucleus cross sections increase with neutrino 

13 



energy and contributions of higher multipole transitions become significant. In general, 
for heavier target nuclei the overall cross sections are more pronounced. In view of the 
modern neutrino detectors based on molybdenium, it is interesting to inspect the results 
of microscopic calculations for the processes induced by supernova neutrinos in the most 
abundant Mo isotopes. In Fig. M the respective cross sections, obtained by folding with 
the Fermi-Dirac distribution (TTTj) for a=0, are shown as a function of temperature for even 
isotopes 92_100 Mo. In accordance with the cross sections shown in Table |TT1 the averaged 
cross sections increase with the number of neutrons in the Mo isotope chain. However, the 
differences between the averaged cross sections are more pronounced at lower temperatures 
due to larger sensitivity of the cross sections to the transitions involved. For example, the 
ratio < ct( 100 Mo) > / < cr( 92 Mo) >=(2.46,1. 31, 1.18,1. 14,1. 13) for the set of temperatures 
T=(2,4,6,8,10) MeV, respectively. 

IV. CONCLUSION 

In summary, modeling of the neutrino-nucleus scattering through the weak neutral cur- 
rent provides important data for simulations of supernova evolution and detector response to 
neutrinos emerging from explosive stellar events. Due to lack of experimental data, it is nec- 
essary to provide independent microscopic insights into the properties of neutrino-induced 
processes, and assess the theoretical uncertainty inherent to implementation of various nu- 
clear effective interactions which determine the transition matrix elements contributing to 
the neutrino-nucleus cross sections in the neutral channel. 

In this work the self-consistent framework for inelastic neutral-current neutrino-nucleus 
scattering is introduced, based on systematic implementation of relativistic nuclear energy 
density functional with density dependent meson-nucleon couplings. The cross sections have 
been formulated using the weak interaction Hamiltonian and nuclear properties of initial and 
excited states are obtained by using the RHB+RQRPA, thus allowing studies of open shell 
target nuclei that necessitate explicit inclusion of the pairing correlations. In order to include 
complete set of natural and unnatural parity excited states, the RQRPA residual interaction 
has been extended using the pion contributions with pseudovector coupling. In the present 
analysis, the neutral-current neutrino-nucleus cross sections have been calculated for the 



set of target nuclei, 12 C, 40 Ar, 56 Fe, and 56 Ni. In addition, in view of the MOON 
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12J and 



NEMO [I4j experiments based on molybdenum detectors, the present study covered the 
respective neutrino-nucleus cross sections in the neutral channel for the most abundant even 
isotopes 92 ~ 100 Mo. In addition to tables |T] and HT1 presented in this work, complete tables of 
all calculated cross sections with small step in neutrino energy are available on request. 

Comparison of the cross sections and their multipole composition appear to be in rea- 
sonable agreement with previous studies, however, some quantitative differences have been 
observed. From the comparison with calculations based on hybrid model and QRPA, the 
present analysis provides an estimate of the theoretical uncertainty in modeling the cross 
sections in the neutral channel due to implementation of various theory frameworks and 
nuclear effective interactions. In the case of 56 Fe, it is shown that the overall cross sections 
exhibit variations, i.e., at some neutrino energies the cross sections based on RQRPA, QRPA 
and hybrid model can differ by a factor «2. This result, together with discrepancies up to 
an order of magnitude between the RQRPA and QRPA calculations shown for 12 C and 40 Ar, 
indicates that in future studies one could critically assess the effect of uncertainties emerg- 
ing from the calculated cross sections on supernova evolution models or neutrino detector 
response. 

The main advantage of the present approach to the neutrino-nucleus cross sections in 
the neutral channel is self-consistent modeling of all relevant transition matrix elements 
involving open-shell nuclei, without any additional adjustments of the model parameters to 
the nuclear target under consideration. In this way, the present study paves the way for 
systematic self-consistent large-scale calculations of stellar neutrino-nucleus scattering in the 
neutral channel. However, this goal would also necessitate further extension of the model, in 
order to include the finite temperature effects in description of nuclei and their excited states 
in the supernova environment. As shown in Refs. 7|, |23j, at finite temperature the cross 
sections become somewhat enhanced at lower neutrino energies. In the forthcoming study, 
the present theoretical framework will be extended to include finite temperature effects, 
typical for the supernova environment. 
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FIG. 1. (Color online) Dependence of the neutrino-nucleus cross sections for the scattering process 
12 C(^ e , fg) 12 C* on the incoming neutrino energy. The cross sections with separate contributions 
from various multipoles J w = — 3 and the total cross sections, including J w = — 5 states, 
are shown. The overall cross sections (stars) are shown in comparison to the QRPA based results 
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FIG. 2. (Color online) The same as in Fig. [TJ but for the scattering process 40 Ar(f e , i/l) 40 Ar*. The 
overall cross sections (stars) are shown in comparison to the QRPA based results (full circles) from 
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FIG. 3. (Color online) The same as in Fig. [H but for the scattering process 56 Fe(^ e , fg) 56 Fe 



,n56tv* 



o 



C\J 



LU 
O 



sum 
+ 


X 56 Ni 


0" 


■■■■■■■■ ^_______-i — — 


1" 

1 + 
-2 + 


— a— __----+'"' a-- 

_^- a— _..»..»(.'"*■■:" 

■ ■■■■■■ ^*r : » , .."J'-* - * 

o ^ — ■• .---"_.-»•;:.-"" .a-- 


. 2" 

3" 


A S^''''\.^ € ^<^ > - - - " " " ■ 

" A " jS s>' „••'• 


"3 + 


■-*•■-/ / ..-'". 




/ ^ 








/ // S *' ..-■- 




/ / / y •* .•■■** 




/ ' •' / m'' _.••••"■" 




/ / ' /' .• .•••" 


_ 


/ ' /'■' •' .-•■'■ 




*' // / ■■•' ..*•'"' 




/ i i' 




I if / .'■ ..-•••■ 




1 ' ft • *■• 


- 


if! / .' 




' if' i .•'•••*" 




.'/ / M' 




/ 1/ ••*'!•'" 




t . S . i . . / . .• i . . . . ' 



10 20 30 40 50 60 70 

E v [MeV] 



80 90 100 



FIG. 4. (Color online) The same as in Fig. [JJ but for the scattering process Ni(z^ e , fg) 56 Ni 
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FIG. 5. (Color online) The same as in Fig. Q3 but for the scattering process 98 Mo(f e , ^) 98 Mo 
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FIG. 8. (Color online) Neutral-current neutrino-nucleus cross sections averaged over the supernova 
neutrino flux in the temperature interval T„ = 2 — 10 MeV. Results are shown for 12 C, 40 Ar, 56 Fe, 
56 Ni and 98 Mo target nuclei. 
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FIG. 9. (Color online) The same as in Fig [8] but for 92,94,96,98,i00 Mo target nuclei 
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